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NOMENCLATURE 


e t - Total energy 

M - Mach number 

P - Static pressure 

q - Magnitude of the velocity 

Re - Inlet reference Reynolds number 

T - Static temperature 

it, v - x, y components of velocity 

a - Absolute flow angle 

/ 3 - Relative flow angle 

p - Density 

r - Shear stress 

D - Rotor rotational speed 

SUBSCRIPTS 

k - Kinetic energy 

rel - Relative reference frame quantity 

t - Stagnation quantity 

T - Turbulent quantity 

w - Wall quantity 

x,y - Paxtial derivative with respect to x, y 
e - Turbulent dissipation 
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ABSTRACT 


Experimental data from jet-engine tests have indicated that unsteady blade-row interaction 
effects can have a significant impact on the efficiency of low-pressure turbine stages. Mea- 
sured turbine efficiencies at takeoff can be as much as two points higher than those at cruise 
conditions. Preliminary studies indicate that Reynolds number effects may contribute to 
the lower efficiencies at cruise conditions. In the current study, numerical experiments have 
been performed to quantify the Reynolds number dependence of unsteady wake/separation 
bubble interaction on the performance of a low-pressure turbine. 

INTRODUCTION 

The need for an improved understanding of efficiency variations in multiple blade row tur- 
bines has motivated engineers to assess the effects of unsteady aerodynamic interaction on 
the performance of turbomachinery components. The two principal types of aerodynamic 
interaction are usually referred to as potential-flow and wake interaction. Potential-flow in- 
teraction results from the variations in the velocity potential or pressure fields associated 
with the blades of a neighboring row, and their effect upon the blades of a given row moving 
at a different rotational speed. This type of interaction is of serious concern when the axial 
spacing between adjacent blade rows is small or the axial flow Mach number is high. Wake 
interaction is the effect upon the flow through a downstream blade row of the vortical and 
entropic wakes shed by one or more upstream rows. 

The focus of the current effort is to use a viscous, unsteady Navier-Stokes analysis to study 
the effects of Reynolds number variation on the interaction of blade wakes and separation 
bubbles in a low-pressure turbine blade passage. The change in the efficiency as a function 
of Reynolds number and separation bubble behavior has been investigated. 


NUMERICAL PROCEDURE 


Governing Equations 

The field equations considered in this investigation are the time-dependent Navier-Stokes 
equations. These equations can be written in non-dimensional form, and in terms of body- 
fitted coordinates, as 


Qr + (F, + Re~ l F v )^ + (Gi + Fe - 1 G„)^ = 0 ( 1 ) 

where the subscripts i refer to inviscid quantities, the subscripts v refer to viscous quanti- 
ties, the subscripts f and rj refer to derivatives with respect to the body-fitted coordinate 
directions, and Re is the Reynolds number. The vector of conserved variables is 


Q = JQ 


( 2 ) 


= J((tQ + £ X Fi + tyG>) 

= j ( VtQ + IxFi + rj y Gi) 


and the inviscid flux vectors are 

Fi 

Gi 
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( 3 ) 

( 4 ) 



The Jacobian of the transformation (from physical Cartesian ( x,y ) to computational curvi- 
linear (£,t?) coordinates) and the other metric quantities are given by 

J = HVn ~ Vi x n ( 5 ) 

£r = J/r ;/ J) £jr = = Vi! Vy ~ ^ (®) 

In the Cartesian coordinate frame, the conserved variable and inviscidflux vectors are defined 
as 

Q = \P,pu,pv,e t \ 

Fi = [pu, pu 2 + P, puv, (e< + P)it] (7) 

Gi = [pu, puv, pv 2 + P, (e t + P)u] 

where p is the density, u and v are the Cartesian velocity components, P is the thermody- 
namic pressure, and e t is the total internal energy. 

The governing equations can be simplified by incorporating the thin-layer assumption [1]. 
The thin-layer assumption assumes that the viscous flux terms parallel to a solid surface are 
negligible with respect to the viscous terms normal to the surface. Thus, for two-dimensional 
turbomachinery applications, the viscous terms in the direction normal to the blade surface 
(77 direction) are retained, 

Qr + (#)< + (Gi + Re~ 1 G v ) n = 0 ( 8 ) 

The viscous flux vector, G v , in Eqn. 8 is given by 


G v = - 


0 

Kiu v + K 2 r) x 
Kiv v + K 2 r) y 

k , [p,-‘(7 - irV), 

+(« 72 )J + k 2 k 3 


( 9 ) 


Ki = fi (q 2 + 7/j, 2 ) 


K 2 = ^ {TJxUr, + TjyV v ) 

( 10 ) 

K 3 = UTfe + VT)y 


q 2 = u 2 + v 2 



where p is the dynamic viscosity, P r is the Prandtl number, and a is the speed of sound. In 
the current study, both the thin-layer and full Navier-Stokes equations have been utilized. 


Solution Algorithm 

The implicit numerical procedure used in this study, consists of a time-marching finite- 
difference scheme [2,3,4]. The procedure is second-order temporally accurate and third-order 
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spatially accurate. The inviscid fluxes are discretized according to either the scheme devel- 
oped by Chakravarthy and Osher [5] or that of Roe [6]. The viscous fluxes axe discretized 
using standard central differences. An alternating direction, approximate-factorization tech- 
nique is used to compute the time-rate changes in the primary variables. Newton sub- 
iterations axe used at each global time step to increase stability and reduce linearization 
errors. For typical unsteady simulations, two Newton sub-iterations are performed at each 
time step. The grid topology for the numerical procedure consists of 0- and H-type overlaid 
grids. 


BOUNDARY CONDITIONS 

Inlet, Exit, and Surface Boundary Conditions 

The inlet, exit, and surface boundary conditions used in this study axe based on a char- 
acteristic analysis of the linearized Euler equations. For viscous flows, it is assumed that the 
computational inlet and exit axe in regions where the flow is predominantly inviscid, and the 
characteristic boundaxy conditions axe retained. No-slip boundary conditions are enforced 
at solid surfaces for viscous flows, along with either a specified wall temperature or heat flux 
distribution. 

In the computational procedure, quasi-two-dimensional characteristic boundaxy condi- 
tions are solved implicitly along with the interior of the computational domain. After each 
time step, fully two-dimensional unsteady boundary conditions can be explicitly applied at 
the inlet and exit boundaries to improve solution accuracy. Further information describing 
the inlet, exit, and surface boundary conditions can be found in Giles [7] and Dorney and 
Davis [3]. 


Zonal Boundary Conditions 

In the numerical analysis, the use of overlaid grids requires the application of zonal 
boundary conditions. Dirichlet conditions, in which the time rate change in the conserved 
vector Q is set to zero, are imposed at the overlaid boundaries of the 0- and H-type grids. 
The flow variables of Q at the zonal boundaries axe explicitly updated after each time 
step by interpolating values from the adjacent grid. Because of the explicit application 
of the zonal boundary conditions, large time steps necessitate the use of more than one 
Newton sub-iteration to maintain time accuracy. The accuracy of the information transfer 
between adjacent grids can also be enhanced by increasing the amount of overlap between 
the 0- and H-grids. The zonal boundaxy conditions axe non-conservative, but the current 
implementation has been shown to yield satisfactory results even for transonic and supersonic 
flows [3, 8]. Further information describing the zonal boundaxy conditions can be found in 
Rai [4, 9]. 


TURBULENCE AND TRANSITION MODELS 

Two turbulence models axe utilized in the numerical procedure; an algebraic model and 
a model based on the solution of two partial differential equations. Several correlations are 
available in the analysis for predicting natural and bubble transition. 
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Baldwin-Lomax Model 

In the algebraic, two-layer Baldwin-Lomax (B-L) turbulence model, the eddy viscosity is 
described by [1] 

_ f pT inner 5 — 5 crossover Hi) 

“ \ pT outer s ^ ^crossover 

where s is the distance normal to the solid surface and s crossover is the smallest value at which 

PT inner PT outer* 

In the inner region, the eddy viscosity is calculated using the Prandtl~Van Driest formu- 

Winner = WV P ( 12 ) 

where l is the mixing length given by 

l = ks 1 — exp(— y + /A + )j (13) 

and k and A + are constants. While formulae exist for modifying k and A + for the presence 
of pressure gradients and surface roughness (see Granville [10] for a review), constant values 
are used here for simplicity. The magnitude of the vorticity, |o>|, can be written in two- 
dimensions, for example, as 


and y + is the law-of-the-wall coordinate 


1 — \ u y v x\ 

(14) 

l_ \j Pw 5 

(15) 

Pw 


calculated using 


K CqpP Fwake FKleb(s') 

(16) 


where K is the Clauser constant, is an additional constant, and F wake is described by 

F wake = min (s ma xF max? Cwk'Smax'J'dif / Fmax) (1"0 

The term F max is the maximum value of F(s) along a given computational grid line normal 
to the surface and . , 

F(s) = s|w| [l - exp (-y + /x4 + )j (18) 

The Klebanoff intermittency factor, F K i e b(-s) is given by 

Fkleb(s) = [l + 5.5((sCKleb)/ 5 max) ] (19) 

and qdii is the difference between the maximum and minimum total velocity in the profile. 
The constants used in the current implementation of the B-L turbulence model are 

A+ = 26 Cop = 1.6 

C K ieb = 0.3 C wk = 0.25 (20) 

k = 0.4 K = .0168 
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Several researchers have developed correlations to modify the constants listed above [11, 12], 
but these relations are usually case-dependent and require significant adjustment by the user. 

One modification that has been made to the B-L model is to apply a smoothing function 
to the predicted turbulent field for separated flow fields. Prudent use of the smoothing 
function has been found to eliminate non-physical gradients in the turbulence field [13]. 

The B-L model has two limitations which are of great consequence in turbomachinery 
simulations. First, the B-L model was designed for steady flow fields and does not contain 
time dependent terms ( i.e ., convection effects cannot be modeled). Thus, the results of 
unsteady simulations using the B-L model must be interpreted with care. Second, the B-L 
model is not dependent on the free stream turbulence level. The free stream turbulence 
level can have a significant impact on the boundary layer and secondary flow development 
in turbomachines. 


k — t Model 


The two-equation k — t model, originally developed by Jones and Launder [14, 15], uses 
transport equations (based on the Navier-Stokes equations) for the kinetic energy (k) and 
turbulent dissipation (e). The k — e equations include the convection and upstream history 
of the turbulence, overcoming the major deficiencies of the B-L model. 

The k — t equations can be written in two dimensions as 


where 


and 


dA dB dC _ J_ (dM_ dN_ ' 
dt dx dy Re \ dx + dy 



B = 


puk 

pue 


pvk 

pve 


M = 

' „ 9fc ' 

N = 

„ dk i 
^ k dy 
Sc 




^Ty J 


S = 


P — peRe 

C^P-^p^Re 


Pk = + Hrfak) 

/*. = (f» + wW 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


P represents the production of kinetic energy and can be written in simplified form as [16] 


P = PT 


(u, 2 + u, 2 ) 


(29) 
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(30) 


The turbulent viscosity is related to the kinetic energy and turbulent dissipation by 

k 2 

V-T = exp- 
and the constants in the model are usually given by [16, 17, 18]: cr fc = 1.0, <r £ = 1.3, 
Ci = 1.35, C 2 = 1.80 and C M = 0.09. The k - e equations are integrated in time using the 
same numerical technique as that used for the Navier-Stokes equations. 

The equations above represent the high Reynolds number form of the k — e equations. The 
equations must be augmented by a near-wall treatment for determining the flow characteris- 
tics in the vicinity of solid surfaces. This is usually accomplished using low Reynolds number 
corrections ( e.g Refs. [15, 19, 20, 21, 22, 23]) or wall functions ( e.g ., Refs. [16, 17, 24, 25]). 
For the simulations presented in this study, the low Reynolds number correction of Chien [20] 
was utilized. 

Although the k — t model includes more relevant physics than the B-L model, it does 
have certain drawbacks which should be noted. First, there are 5 constants which must be 
specified (more if low Reynolds number corrections are used). Second, boundary conditions 
must be specified at the inlet, exit and solid walls. Perhaps the most troublesome boundary 
conditions are at the inlet, where the free stream turbulence level and a length scale are 
normally specified. While reasonable estimates can be made for free stream turbulence level, 
the appropriate length scale is difficult to measure experimentally and requires significant 
intuition on the part of the user [26, 27]. 

Transition Models 

There has been a considerable amount of interest over the past decade in developing reli- 
able transition models which provide accurate results in zero, favorable and adverse pressure 
gradients. Several research groups are actively involved in developing transition models for 
turbomachinery applications based on a combination of analytical, experimental and compu- 
tational work [28, 29, 30, 31, 32, 33]. Although these newer correlations show great promise, 
the majority of industrial codes still incorporate simpler transition models. 

The current numerical procedure contains four different models for predicting the lo- 
cation of natural transition. The four models represent correlations of data from flat plate 
experiments, and therefore yield varying results depending on the curvature of the blades. In 
all the models, the onset of transition is calculated by determining when the local Reynolds 
number (based on the momentum thickness) exceeds a specified value. 


Hall and Gibbings Model 

The momentum thickness Reynolds number signifying the start of transition in the model 
developed by Hall and Gibbings [34] is given by 

Re e = 190 + exp (6.88 - 1.03r) (31) 

where r is the free stream turbulence level measured as a percentage of the free stream 
velocity. This transition model is valid for zero pressure gradients and turbulence levels up 
to approximately 9 percent. 
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Seyb and Singh Model 


The momentum thickness Reynolds number indicating the beginning of transition based 
on the model of Seyb [35] and Singh [36] is determined by 


Reg = 


1000 


+ 10 


1.2 + 0.7r 

where \$ is the pressure gradient parameter defined by 

e 2 du - 


( 0.09 + A* \ 

\ 0.0106 + 0.036r ) 


2.62 


(32) 


Afl = 


v ds 


(33) 


In Eqn. 33, 6 is the boundary layer momentum thickness, v is the kinematic viscosity, [/«> is 
the velocity at the edge of the boundary layer, and s is a distance parameter measured from 
the leading edge stagnation point. 


Dunham Model 

The momentum thickness Reynolds number corresponding to the start of transition in 
the model developed by Dunham [37] is 

Re e = [0.27 + 0.73 exp (-0.8r)] (550 + 68 ° 0i . ) (34) 

The Dunham model has been widely used, including the study presented in Ref. [38]. 


Abu-Ghannam and Shaw Model 


The momentum thickness Reynolds number indicating the beginning of transition based 
on the work of Abu-Ghannam and Shaw [39] is given by 


where 


Reg = 163 + exp 


F(X,) - 


F(M r 

6.91 


F(A») = 6.91 + 12.75A, + 63.64(A») ! A, < 0 
F(A # ) = 6.91 + 2.48 A* - 12.27(A,) ! A, > 0 


(35) 


(36) 

(37) 


The Abu-Ghannam and Shaw model has been successfully used in several turbomachinery 
studies, including Refs. [40] and [41]. Note, all the transition simulations in this study were 
performed using the Abu-Ghannam and Shaw model. 

For all the transition models, the end of the transition region (i.e., fully turbulent flow) 
is identified by the position where the local Reynolds number, Rez, exceeds the following 
value [39] 

Re L = 31.8(Re e ) 16 (38) 
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where the momentum thickness Reynolds number is equal to its value at the start of transi- 
tion. In the region between the beginning and end of transition the intermittency factor, <r, 
is varied according to the theory of Dwahan and Narasimha [42] 

<7 = 0 x < X tr b (39) 

a = 1 — exp (— 4.64£ 2 ) x tr i , < x < xt re (40) 

a — 1 X > Xfre (41) 

where x tr b denotes the beginning of transition, Xt re denotes the end of transition, and 


X Xf r b 
Xfre % trb 


(42) 


In cases where separation bubble separation occurs, the correlation developed by Roberts 
is implemented [43, 44]. All the transition models described in this section were developed 
from two-dimensional data and are therefore suspect in three-dimensional flow simulations; 
there are currently no reliable correlations for general three-dimensional turbomachinery flow 
fields. 


NUMERICAL SIMULATIONS 

The geometry studied in this investigation is the Pratt & Whitney “PAK B” geometry. 
The characteristics of this cascade include: 

• Ti— 530°R, Pi = 14.7psia 

• Mi = 0.0897 (based on a Q-1D analysis), M 2 = 0.15 

• c*i = 35°, <*2 = —60° (measured from the axial direction) 

• Pi/Pti ~ .9844 (based on a Q-1D analysis) 

• Re = 27, 000 - 200, 000 based on the chord and the inlet flow variables 

• Pitch-to-chord ratio = .88558 

• Adiabatic walls 

Laminar, turbulent, and transitional simulations have been performed to better understand 
the role of low Reynolds number effects on unsteady separation bubbles. 

Baseline PAK B Cascade Simulations 

Two computational grids were used in the cascade simulations; the coarser grid contained 
251 x 51 (streamwise x tangential) points in the O-grid and 180 x 45 points in the H-grid for 
a total of 20,901 grid points. The average value of y + , the non-dimensional distance of the 
first grid point above the surface, was approximately equal to 0.17. The finer computational 
grid was generated with 281 x 51 points in the O-grid and 240 x 45 in the H-grid for a total 
of 25,131 computational grid points.The average value of y + was again about y + = 0.17. 
A comparison of the coarse and fine computational grids is presented Table 1. Typical 
simulations require 0.00029 seconds of CPU time per grid point per iteration on a DEC 
3000-300 workstation. 
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Re = 40, 000, B-L, Turbulent Flow 


Figure 2 illustrates the surface pressure distributions predicted using the coarse and 
fine grids. The two solutions show fair agreement, with small differences noticeable on 
the pressure surface and upstream of the throat on the suction surface. Figure 3 contains 
the surface isentropic Mach number distributions from the coarse and fine-grid simulations. 
The design Mach number distribution is also included in Fig. 3. The results of the fine-grid 
simulation show very good agreement with the design data. The coaxse grid solution exhibits 
very good agreement with the design data on the suction surface of the blade, but slightly 
underpredicts the Mach number on the pressure surface upstream of approximately 65% of 
the axial chord. The predicted turbulent skin friction distributions are shown in Fig. 4. 
The fine-grid skin friction distribution is nearly identical to the coarse-grid distribution, 
including the suction surface separation location. The momentum thickness distributions 
from the two simulations were also similar. Unless otherwise noted, subsequent figures will 
contain fine-grid results. 

Table 2 contains a comparison of the average inlet and exit quantities from the coarse 
and fine-grid simulations. Note, the low-speed nature of the flow causes minor changes in 
the inlet or exit Mach number/stagnation quanitities to have a significant impact on the 
predicted losses. 


Re = 40, 000, Laminar Flow 

Both the coarse and fine-grid laminar flow solutions exhibited a large, unsteady suction- 
surface separation region. The boundary layer velocity profiles were investigated at different 
locations to ascertain that an adequate number of grid points were contained in the boundary 
layer. Figure 5 illustrates velocity profiles at 5% and 65% of the axial chord along the suction 
surface. Near the leading edge the boundary layer region contains 16 grid points, while near 
the separation point the boundary layer contains 25 grid points. Figure 6 shows velocity 
profiles at 5% and 75% of the axial chord along the pressure surface. Similar to the suction 
surface, the boundary layer near the leading edge contains 16 grid points and the boundary 
layer near the trailing edge contains approximately 25 grid points. The resolution of the 
boundary layer appears to be adequate over the entire blade surface. Therefore, the suction 
surface unsteady flow separation is probably not caused by an insufficient number of grid 
points. 

Figures 7 and 8 illustrate instantaneous Mach number and entropy contours, respectively, 
from the fine-grid simulation. Figures 7 and 8 are similar to the results of the coarse-grid 
simulation. Figure 9 and 10 show the pressure and skin friction envelopes, respectively. The 
amplitude of the unsteadiness along the suction surface (see Fig. 9) and the time-averaged 
separation location (see Fig. 10) were nearly identical to the coarse-grid results. 

Table 2 contains the time-averaged inlet and exit flow quantities from the coarse and 
fine-grid laminar flow simulations. 


Re = 40,000, B-L, Transitional Flow 

Coarse and fine-grid simulations were performed with the onset of transition prescribed to 
be at 88% of the axial chord on the suction surface. The transition region extended over 5% of 
the axial chord. The beginning and end of transition were determined by time-averaging the 
unsteady transition locations predicted using the Abu-Ghannam/Shaw transition model [39]. 
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Figures 11 and 12 contain instantaneous Mach number and entropy contours, respectively, 
from the fine-grid simulation. Fixing the transition location eliminates the unsteadiness 
observed when the position is allowed to vary with time. In addition, the convection of 
the airfoil wake appears to be more accurately resolved using the refined grid (see Fig. 12). 
The time-averaged surface isentropic Mach number distributions predicted in the coarse and 
fine-grid simulations are shown in Fig. 13. The design Mach number distribution is also 
included in the figure. The coarse and fine-grid solutions are similar and exhibit excellent 
agreement with the experimental data. Figure 14 illustrates the skin-friction envelopes for 
the transition flow conditions. Again, the coarse and fine-grid solutions show very good 
agreement. 


Re = 120,000, B-L, Transitional Flow 

Figure 15 illustrates the predicted and design distributions of the surface (isentropic) 
Mach number for Re = 120,000. The predicted results in Fig. 15 were obtained using 
the average values of the transition locations predicted by the Abu-Ghannam and Shaw 
model; the suction surface transition location was set at x* r = 0.88 and the pressure surface 
transition location was specified to be at x tr = 0.25. The predicted and design Mach number 
distributions exhibit close agreement. Figure 16 contains entropy contours for the simulation 
with fixed transition. The flow field does not display the unsteady shedding from the suction 
surface displayed at Re = 40, 000. In addition, the wakes remain undisturbed as they are 
convected to the exit boundary. The baseline inlet and exit flow quantities, as well as the 
cascade loss, are presented in Table 3. 


Re = 200,000, B-L, Transitional Flow 

A history of the suction surface transition locations predicted using the Abu-Ghannam 
and Shaw transition model is shown in Fig. 17. The average transition location is at approx- 
imately 80% of the axial chord (x tr = 0.80). Thus, the location of transition on the suction 
surface is 8% upstream of where it was predicted at Re = 120,000. The variations in the 
location of transition shown in Fig. 17 appear to be time-periodic. Since no unsteady shed- 
ding was observed in the flow field, the Abu-Ghannam and Shaw transition model may be 
getting trapped in a limit cycle. The transition location on the pressure surface of the blade 
was predicted to be at the leading edge (i.e., the flow on the pressure surface is completely 
turbulent). A comparison of the predicted (using fixed transition) and design surface Mach 
number distributions are illustrated in Fig. 18. The two sets of data show good agreement, 
with only small discrepancies visible on the pressure surface near 25% of the axial chord. 
Similar to the Re = 120, 000 results, the entropy contours (see Fig. 19) do not indicate 
unsteady shedding from the suction surface of the blade. The baseline inlet and exit flow 
quantities, as well as the cascade loss, are presented in Table 3. 

Effects of Turbulence Modeling Parameters - k - e 

A series of simulations have been performed for the PAK B isolated cascade geometry 
using the k — e turbulence model and a free stream turbulence level of Tu = 3.0%. The 
Reynolds number was set at Re = 80, 000. The simulations utilized the turbulence field 
predicted using the Baldwin-Lomax model as the initial solution. Figure 20 contains the 
skin friction distributions predicted using the B-L turbulence model [47], and the k - e 
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model with different values of the free stream turbulent viscosity. Setting the free stream 
turbulent viscosity equal to unity creates too much dissipation and the skin friction behaves 
in a laminar manner, as evidenced by a large separated flow region on the suction surface 
of the airfoil. Increasing the free stream turbulent viscosity to a non-dimensional value of 
10.0 decreases the dissipation, and the skin friction distribution shows good agreement with 
the results obtained using the B-L model. The results suggest that the predicted results are 
strongly dependent on the free stream value of the turbulent viscosity (and consequently e). 
Figure 21 displays a comparison of the momentum thickness distributions predicted using 
the two turbulence models. The results obtained using the k — e model with /z too = 10.0 show 
fair agreement with the results obtained using the B-L model. The k — e model produces 
a greater momentum thickness on the pressure surface, while the B-L model produces a 
greater momentum thickness on the suction surface. Figures 22-24 show the Mach number, 
entropy and turbulent viscosity contours, respectively, predicted using the k — e model with 
Ht<x> = 10.0. The contours are much smoother then those predicted using the B-L model, 
especially those for the turbulent viscosity. Figure 25 illustrates the predicted and design 
distributions of the surface (isentropic) Mach number. The predicted results and the design 
data exhibit good agreement. 

The simulations with using //Too = 10-0 were repeated with the free stream turbulence 
level increased to Tu = 6.0%. Figures 26 and 27 illustrate comparisons of the skin friction and 
momentum thickness distributions, respectively, with free stream turbulence levels of 3.0% 
and 6.0%. The skin friction distributions from the two simulations are similar (see Fig. 26), 
except that Tu = 6.0% produces slightly smaller values over the first 50% of the axial chord 
on the suction surface. The momentum thickness distributions are comparable on the suction 
surface of the airfoil (see Fig. 27). In the cove region on the pressure surface, however, the 
lower free stream turbulence level generates a much larger momentum thickness. Figures 28- 
30 show the Mach number, entropy and turbulent viscosity contours, respectively. The three 
figures are similar to those obtained for Tu = 3.0% (see Figs. 22-24), which suggests that 
the solution is more dependent on the initial free stream value of e than on the free stream 
turbulence level. 


PAK B Cascade Wake-Passing Simulations 

Re = 40, 000, LSRR Stator Wake, B-L, Transitional Flow 

To accurately model the flow physics of a turbine stage, a wake definition was obtained 
from a numerical simulation for the first-stage stator of the Large Scale Rotating Rig (LSRR) 
turbine [48]. In the unsteady numerical simulations utilizing the LSRR stator wake, the 
perturbations in the velocity components, pressure and density have been included in the 
definition of the unsteady inlet boundary conditions. 

The LSRR stator geometry was rotated 32 degrees for use in the current investigation (see 
Fig. 31), such that the flow would be turned approximately 135 degrees by the combination 
of the stator and the PAK B turbine blade. The wake velocity profile generated by the 
numerical simulation for the stator is shown in Fig. 32. The wake width is approximately 
25% of the stator gap, and the wake deficit is approximately 10% of the free stream velocity. 
The unsteady simulations were performed assuming one upstream stator wake for every 
two rotor blade passages (i.e., interblade phase angle of a = 7 r), and the reduced temporal 
frequency was set to k = 5. One unsteady simulation was performed with a fixed value for 
the transition location on each surface of the blade, and one simulation was performed using 
the Abu-Ghannam and Shaw model to predict the location of natural transition [39] and 
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Roberts correlation to predict bubble transition [43, 44]. 

Figure 33 illustrates perturbation (i.e., instantaneous minus time-averaged) vorticity con- 
tours for the simulation with a fixed location for transition. The vorticity associated with 
the stator wake is generally smaller than the values in the viscous layer near the blade, but 
it can still be used to highlight the path of the wake upstream of the blade row. Because the 
pressure and density are allowed to vary through the wake, the entropy can be used to track 
the movement of the wakes through the entire blade passage. Figures 34 and 35 illustrate 
instantaneous entropy contours from the simulations with fixed and predicted (floating) tran- 
sition, respectively. The stretching of the wake between the pressure and suction surfaces 
of adjacent blades is evident in the blade passages. The interaction between the stator and 
blade wakes downstream of the blade passages is also discernible. Comparing Figs. 34 and 
35, it is observed that the simulation with floating transition exhibits more unsteadiness in, 
and downstream of, the blade passage. 

The predicted skin friction envelopes from the two unsteady simulations are shown in 
Figs. 36 and 37, respectively. The results using floating transition (see Fig. 37) exhibit more 
unsteadiness along the suction surface of the blade than the results with fixed transition. 
Moreover, the results for fixed transition indicate that the flow downstream of the time- 
averaged separation point intermittently reattaches, while the results for floating transition 
do not show reattachment during any portion of the wake-passing cycle. Along the majority 
of the pressure surface the two solutions show good agreement. Differences are noticeable 
near 25% of the chord, where the case with floating transition appears to be closer to sep- 
aration, and near the trailing edge where the case with floating transition displays more 
unsteadiness. 

Figures 38 and 39 contain the momentum thickness envelopes from the simulations with 
fixed and floating transition, respectively. The predicted results for the pressure surface are 
s imilar in both simulations, except that the simulation with floating transition shows smaller 
values near 25% of the chord (corresponding the differences in the skin friction distributions 
noted above). Along the suction surface of the blade the simulation with floating transition 
shows more unsteadiness, especially near the trailing edge. 

The variation in the flow angle exiting the blade row is illustrated in Fig. 40. The case 
with fixed transition produced a periodic variation of approximately 0.35 degrees in the exit 
flow angle. The simulation with floating transition generated an exit flow angle variation of 
almost 1 degree, and was only marginally periodic. 

In the simulation with floating transition, the temporal variation of the transition loca- 
tions were recorded (see Figs. 41 and 42). In both figures, a value of 0 corresponds to the 
leading edge and a value of 1.0 corresponds to the trailing edge. Along the suction surface 
of the blade (see Fig. 41), the transition point is located at, or near, the trailing edge for 
the majority of the wake passing cycle. During portions of the cycle, the transition location 
moves forward to approximately 75% of the axial chord. When transition was predicted 
upstream of the trailing edge, it was due to bubble transition. Along the pressure surface of 
the blade (see Fig. 42), the location of transition varied between 25% of the chord and the 
trailing edge. Similar to the suction surface results, bubble transition was observed at all 
times along the pressure surface of the blade. When transition was predicted near 25% of the 
chord, it was due to cove separation. An interesting phenomenon was observed when tran- 
sition was predicted upstream of the trailing edge. For several time steps bubble transition 
would be predicted at a given location. During this time, the eddy viscosity added to the 
flow reduced or eliminated the separation bubble. After the separation bubble disappeared 
the transition model did not detect transition until the trailing edge, allowing the separation 
bubble to develop once more. This cycle was repeated several times (see Fig. 43). A similar 
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process was also observed on the suction surface of the blade. 

Table 4 contains the average inlet and exit flow quantities from the baseline simulation, 
and the two simulations utilizing the LSRR stator wake. The time-averaged Mach numbers 
and flow angles show very good agreement with the baseline results. The time-averaged losses 
from both unsteady simulations are greater than the baseline loss. In addition, the losses 
predicted using floating transition are greater than those predicted using fixed transition. 


Re = 120, 000, LSRR Stator Wake, B-L, Transitional Flow 

Figures 44 and 45 illustrate instantaneous entropy contours from the simulations with 
fixed and predicted (floating) transition, respectively. The stretching of the wakes between 
the pressure and suction surfaces of adjacent blades is evident in the blade passages. Down- 
stream of the blade passage there are two distinct wakes, the one specified at the compu- 
tational inlet and the one generated by the PAK B blade row. Comparing Figs. 44 and 
45 it is observed that the two solutions are similar, except that the simulation with fixed 
transition exhibits somewhat more unsteadiness near the blade trailing edge. Note, in the 
simulation using the Abu-Ghannam and Shaw model the location of the suction surface 
transition point moved from 88% of the chord to the trailing edge, while the pressure surface 
transition location moved from 25% of the chord to the leading edge. 

The predicted unsteady skin friction envelope for the simulation with a fixed transition 
location is shown in Fig. 46. The skin friction exhibits a moderate amount of unsteadiness 
along the suction surface of the blade up to approximately 88% of the chord, where tran- 
sition occurs. Downstream of transition a larger amount of unsteadiness is observed. In a 
time-averaged sense the flow along the suction surface of the blade remains attached, but 
intermittent flow separation is detected from approximately 73% of the axial chord (just 
downstream of the throat) to the trailing edge. Along the pressure surface of the blade, the 
flow contains cove separation from approximately 10% to 35% of the chord. The skin friction 
displays little unsteadiness along the pressure surface, except in the separated flow region. 

Figure 47 contains the skin friction envelope for the simulation in which transition was 
predicted. The skin friction shows a considerable amount of unsteadiness over the entire suc- 
tion surface of the blade. Unlike the simulation with fixed transition, there is no intermittent 
flow separation on the suction surface. The time-averaged flow on the pressure surface of 
the blade remains attached, but intermittent separation is again observed in the cove region. 
The differences observed in the skin friction distributions illustrated in Figs. 46 and 47 are 
probably due to the movement of the suction and pressure surface transition points in the 
simulation in which transition was predicted. 

Figure 48 contains the unsteady momentum thickness envelope from the simulation with 
a fixed transition location. There is very little unsteadiness on the suction surface upstream 
of the intermittent separation point. Much larger values of the momentum thickness are 
present from the approximately 75% of the axial chord to the trailing edge. The momentum 
thickness shows significant excursions around the time-averaged values along the pressure 
surface of the blade. The largest values on the pressure surface correspond to the region of 
cove separation discussed above. 

The unsteady momentum thickness envelope for the simulation with floating transition 
is shown in Fig. 49. The unsteadiness about the time-averaged values gradually increases 
with distance from the leading edge along the suction surface of the airfoil. The minimum, 
maximum, and time-averaged values are all greater than the corresponding values for a fixed 
transition location because the flow has become laminar over the entire suction surface. 
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Large levels of unsteadiness are observed on the pressure surface, especially in the cove 
region. Overall, however, the pressure surface values are smaller than the corresponding 
values for fixed transition. 

The temporal variation in the exit flow angle is illustrated in Fig. 50. Both flow angle 
distributions exhibit some high-frequency unsteadiness (perhaps due to trailing-edge vortex 
shedding and/or the movement of the transition location), in addition to unsteadiness at 
the wake-passing frequency. The average exit flow angle from the simulation with floating 
transition is approximately 0.5 degrees greater than that for fixed transition. 

Table 5 contains the average inlet and exit flow quantities from the baseline simula- 
tion (fixed transition), and the two simulations utilizing the LSRR stator wake. The time- 
averaged Mach numbers and flow angles from the unsteady simulation with fixed transition 
show very good agreement with the baseline results. The time-averaged loss from the un- 
steady simulation with fixed transition, however, is over 20% less than in the baseline sim- 
ulation. The unsteady simulation with floating transition displays less turning and higher 
losses than either of the other two simulations. 


Re = 200,000, LSRR Stator Wake, B-L, Transitional Flow 

Figures 51 and 52 illustrate instantaneous entropy contours from the simulations with 
fixed and predicted (floating) transition, respectively. The two solutions compare favorably 
throughout the blade row. Note, in the simulation using the Abu-Ghannam and Shaw 
model the location of the suction surface transition point moved from 80% of the chord to 
the trailing edge, while the pressure surface transition location remained at the leading edge. 

The predicted unsteady skin friction envelope for the simulation with a fixed transition 
location is shown in Fig. 53. Also included in Fig. 53 is the baseline skin friction distribution. 
The skin friction remains relatively constant along the suction surface of the blade up to 
approximately 65% of the chord, where it rapidly increases. This location corresponds to 
the throat of the blade passage. From approximately 70% of the chord to the trailing edge 
the skin friction gradually decreases. The baseline and time-averaged distributions are nearly 
identical over the entire suction surface of the blade, although the nature of the distributions 
differs from that seen at lower Reynolds numbers. No flow separation occurs along the suction 
surface of the blade. Along the pressure surface of the blade the skin friction approaches 
zero, but the flow does not separate. The skin friction displays little unsteadiness along the 
pressure surface of the airfoil. 

Figure 54 contains the skin friction envelope for the simulation in which transition was 
predicted, along with the corresponding baseline results. The skin friction shows a con- 
siderable amount of unsteadiness over the entire suction surface of the blade. In addition, 
the baseline and time-averaged solutions are significantly different, which may be caused by 
the movement of transition to the trailing edge in the unsteady simulation. The baseline 
distributions displayed in Figs. 53 and 54 show good agreement up to the throat, but are 
different downstream of the throat. Allowing transition to float in the baseline simulation 
caused the generation of a separation bubble just downstream of the throat. No separation 
is detected on either surface of the airfoil in the unsteady simulation. 

Figures 55 and 56 contain the unsteady momentum thickness envelopes for the simula- 
tions with fixed and floating transition, respectively. In Fig. 55 the suction surface values 
remain relatively constant, with little unsteadiness, up to the throat of the passage. Down- 
stream of the throat the momentum thickness increases, as do the excursions around the 
time-averaged values. There is significant unsteadiness in the momentum thickness on the 
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pressure surface, especially in the cove region of the airfoil. In general, the momentum thick- 
ness values shown in Fig. 55 axe smaller than those for Re = 120, 000 (see Fig. 48). In Fig. 56, 
the suction surface distribution exhibits more unsteadiness than with fixed transition (see 
Fig. 55). The time-averaged values of the momentum thickness on the suction surface are 
greater than those shown in Fig. 55, which is expected because the transition point has been 
shifted to the trailing edge. The pressure surface momentum thickness distribution agrees 
closely with that for fixed transition because in both cases the flow is fully turbulent. 

A time history of the exit flow angle is illustrated in Fig. 57. The two distributions show 
very good agreement, although there is somewhat more high-frequency unsteadiness in the 
simulation with fixed transition. 

Table 6 contains the average inlet and exit flow quantities from the baseline simula- 
tion (fixed transition), and the two simulations utilizing the LSRR stator wake. The time- 
averaged Mach numbers and flow angles from both unsteady simulations show very good 
agreement with the baseline results. The time-averaged loss from the unsteady simulation 
with fixed transition is approximately 8% less than in the baseline simulation, while the loss 
in the unsteady simulation with floating transition is approximately 13% less than in the 
baseline simulation. 


PAK B Turbine Stage Simulations - Re = 27, 000- Re = 200,000 

A turbine stage was produced using the PAK B geometry for both the vane and rotor 
airfoils. A l-vane/2-rotor blade count ratio was assumed, and the axial gap between the vane 
and blade was specified to be 30% of the rotor chord. Figure 58 shows the computational 
grid topology, which contains 75,393 grid points. The rotational speed of the rotor was 
specified such that the time-averaged relative frame flow angle entering the rotor passage 
was approximately equal to that used in the isolated blade row simulations (fS 2 = 35°). 

A matrix of numerical simulations have been performed for the PAK B turbine stage 
configuration described above. Turbulent, floating transition and fixed transition simulations 
have been performed for Reynolds numbers of Re = 27,000, Re = 40,000, Re = 60,000, 
Re = 80,000, Re = 120,000 and Re = 200,000 based on the inlet flow conditions. These 
Reynolds numbers correspond to a range of Re = 40, 000 and Re — 300, 000 based on the exit 
flow conditions. The simulations were performed using both the B-L and k — e turbulence 
models. To facilitate interpretation of the results, Tables 7-12 contain the time-averaged flow 
quantities at the inlet and exit of each blade row for the turbulent and transitional simulations 
at each Reynolds number. All solutions were time-averaged over two blade-passing cycles, 
although the extent of flow separation and the tabulated results suggest that the solution 
may need to be averaged over additional blade-passing cycles. Note, the efficiencies show 
considerable scatter; this is due in part to the low Mach number flow (small changes in P t 
or T t equates to large changes in the efficiency) and the need for more blade-passing cycles 
in the time-averaging process. 

Figures 59-61 illustrates the vane and rotor flow turning angles for fully turbulent, float- 
ing transition and fixed transition flow assumptions (using the Baldwin- Lomax turbulence 
model) as a function of the Reynolds number. Figure 62 contains the variation of the vane 
and rotor suction-side transition locations as a function of the Reynolds number. The turning 
angles and transition locations are directly related to the size of the suction-surface sepa- 
ration bubbles. The transition from laminar to turbulent flow always occurred as bubble 
transition, not natural transition. In addition, for all the cases tested the time-averaged 
suction-surface transition location was close to the trailing edge (indicating laminar flow 
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over the majority of the suction surface). The skin friction envelopes corresponding to the 
turbulent and transitional solutions at the five Reynolds numbers are shown in Figs. 63-84, 
respectively. The momentum thickness distributions for the different test conditions are 
shown in Figs. 85-106, respectively. 

Samples of instantaneous entropy and Mach number contours from turbulent flow simu- 
lations at Re = 40,000 and Re = 120,000 are shown in Figs. 107-110. 


CONCLUSIONS 

A numerical study of the losses associated with low Reynolds number flow has been 
conducted for a turbine cascade and turbine stage. Laminar, transitional and turbulent flow 
simulations were performed (over a range of Reynolds numbers) using algebraic and two- 
equation turbulence models, as well as algebraic transition models. The predicted results 
indicate increased losses as the Reynolds number is decreased, and highlight the need for 
improved turbulence and transition modeling. In addition, at the lower Reynolds numbers 
(less than approximately Re = 60, 000) 

• the flow field is inherently unsteady 

• the algebraic turbulence model begins to break down and requires significant user 
intervention 

• the algebraic transition models can get trapped in limit cycles 
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Table 1: Coarse and fine grid topologies 
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Table 2: Average inlet and exit flow quantities for PAK B blade at Re = 40, 000 
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Table 3: Inlet and exit flow quantities for the PAK B blade at Re = 120, 000; 200, 000 
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Table 4: Inlet and exit flow quantities for the PAK B blade at Re = 40, 000, transitional, 
LSRR stator wake 
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Table 5: Inlet and exit flow quantities for the PAK B turbine blade at Re = 120, 000 
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Table 6: Inlet and exit flow quantities for the PAK B turbine blade at Re = 200, 000 
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Table 7: Inlet and exit flow quantities for the PAK B turbine stage at Re = 27, 000 
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Table 8: Inlet and exit flow quantities for the PAK B turbine stage at Re = 40, 000 



Table 9: Inlet and exit flow quantities for the PAK B turbine stage at Re = 60, 000 
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Table 10: Inlet and exit flow quantities for the PAK B turbine stage at Re = 80, 000 
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Table 11: Inlet and exit flow quantities for the PAK B turbine stage at Re = 120,000 
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Table 12: Inlet and exit flow quantities for the PAK B turbine stage at Re = 200, 000 
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Figure 1: Computational grid for the PAK B cascade 
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Figure 2: Surface pressure distribution for PAK B turbine - Re = 40, 000, turbulent 



Figure 3: Surface (isentropic) Mach number distribution for PAK B turbine - Re = 40,000, 
turbulent 
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Figure 6: Pressure surface velocity profiles for PAK B turbine - Re = 40, 000, laminar 
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Figure 7: Instantaneous Mach number contours for PAK B turbine - Re = 40,000, laminar 
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Figure 8: Instantaneous entropy contours for PAK B turbine - Re = 40, 000, laminar 
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Figure 10: Skin-friction envelope for PAK B turbine -Re — 40, 000, laminar 
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Figure 11: Instantaneous Mach number contours 
tional 
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Figure 12: Instantaneous entropy contours for PAK B turbine - Re = 40,000, transitional 
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Figure 13: Surface (isentropic) Mach number distribution for PAK B turbine - Re = 40, 000, 
transitional 



Figure 14: Skin friction distribution for PAK B turbine - Re = 40,000, transitional 
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Figure 15: Surface (isentropic) Mach number distribution for the PAK B 
Re = 120, 000, fixed transition 
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Figure 16: Instantaneous entropy contours for the PAK B turbine blade - Re — 120,000, 
fixed transition 
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Figure 17: Suction surface transition location for the PAK B turbine blade - Re = 200,000 
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Figure 18: Surface (isentropic) Mach number distribution for the PAK B turbine blade - 
Re = 200,000, fixed transition 
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Figure 19: Instantaneous entropy contours for the PAK B turbine blade - Re = 200,000, 
fixed transition 
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Figure 20: 
80,000; Tu 


Figure 21: 
blade - Re 



Turbulence model effects on the skin friction for PAK B turbine blade - Re = 
= 3% 



Turbulence model effects on the momentum thickness for the PAK B turbine 
= 80, 000; Tu = 3% 
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Figure 23: Entropy contours for the PAK B turbii 
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Figure 25: Surface (isentropic) Mach number - Re = 80,000; Tu = 3%; UToo — 10.0 



Figure 26: Free stream turbulence effects on the skin friction for PAK B turbine blade - 
Re = 80,000; Tu =3%, 6% 
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Figure 27: Free stream turbulence effects on the momentum thickness for the PAK B turbine 
blade - Re = 80,000; Tu =3%, 6% 
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Figure 28: Mach number contours for the PAK 
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B turbine blade - Re = 80,000; Tu = 6% 
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Figure 29: Entropy contours for the PAK B turbine blade - Re = 80,000; Tu = 
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Figure 30: Turbulent viscosity contours for the PAK B turbine blade - Re = 80,000; Tu 

6%; /it oo = 10.0 
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Figure 31: Modified LSRR geometry used to generate wake 
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Figure 33: Perturbation vorticity contours for the PAK B turbine - Re = 40,000, LSRR 
stator wake, fixed transition 
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Figure 34: Instantaneous entropy contours for the PAK B turbine - Re = 40, 000, LSRR 
stator wake, fixed transition 
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Figure 35: Instantaneous entropy contours for the PAK B turbine - Re = 40,000, 
stator wake, floating transition 
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Figure 36: Unsteady skin friction envelope for the PAK B turbine - Re = 40, 000, LSRR 
stator wake, fixed transition 
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Figure 37: Unsteady skin friction envelope for the PAK B turbine - Re = 40,000, LSRR 
stator wake, floating transition 



Figure 38: Unsteady momentum thickness envelope for the PAK B turbine - Re = 40, 000, 
LSRR stator wake, fixed transition 
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Figure 39: Unsteady momentum thickness envelope for the PAK B turbine -Re — 40,000, 
LSRR stator wake, floating transition 
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Figure 40: Exit flow angle variation for the PAK B turbine - Re = 40, 000, LSRR stator 
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Figure 41: Time- variation of suction surface transition location - Re = 40,000, LSRR stator 
wake 
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Figure 42: Time- variation of pressure surface transition location - Re = 40, 000, LSRR stator 
wake 
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Figure 44: Instantaneous entropy contours for the PAK B turbine blade - Re = 120,000, 
LSRR stator wake, fixed transition 
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Figure 45: Instantaneous entropy contours for the PAK B turbine blade - Re = 120, 000, 
LSRR stator wake, floating transition 
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Figure 46: Unsteady skin friction envelope for the PAK B turbine blade - Re = 120, 000, 
LSRR stator wake, fixed transition 



Figure 47: Unsteady skin friction envelope for the PAK B turbine blade - Re = 120,000, 
LSRR stator wake, floating transition 
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Figure 48: Unsteady momentum thickness envelope for the PAK B turbine blade - Re 
120, 000, LSRR stator wake, fixed transition 



Figure 49: Unsteady momentum thickness envelope for the PAK B turbine blade - Re 
120, 000, LSRR stator wake, floating transition 
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Figure 50: Exit flow angle variation for the PAK B turbine blade - Re = 120, 000, LSRR 
stator wake 
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Figure 51: Instantaneous entropy contours 
LSRR stator wake, fixed transition 
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Figure 52: Instantaneous entropy contours for the PAK B turbine blade - Re 
LSRR stator wake, floating transition 
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Figure 53: Unsteady skin friction envelope for the PAK B turbine blade - Re = 200, 000, 
LSRR stator wake, fixed transition 



Figure 54: Unsteady skin friction envelope for the PAK B turbine blade - Re = 200,000, 
LSRR stator wake, floating transition 
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Figure 55: Unsteady momentum thickness envelope for the PAK B turbine blade - Re 
200, 000, LSRR stator wake, fixed transition 



Figure 56: Unsteady momentum thickness envelope for the PAK B turbine blade - Re 
200,000, LSRR stator wake, floating transition 
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Figure 57: Exit flow angle variation for the PAK B turbine blade - Re = 200,000, LSRR 
stator wake 




Figure 58: Computational grid for the PAK B turbine stage 
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Figure 59: Flow turning angles; Baldwin-Lomax - turbulent 
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Figure 60: Flow turning angles; Baldwin-Lomax - floating transition 
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Figure 61: Flow turning angles; Baldwin- Lomax - fixed transition 



Figure 62: Suction surface transition locations 
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Figure 63: Skin friction envelope for the PAK B turbine stage - Re — 27 , 000; B-L, turbulent 



Figure 64: Skin friction envelope for the PAK B turbine stage - Re — 27,000; B-L, floating 
transition 
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Figure 65: Skin friction envelope for the PAK B turbine stage - Re = 27,000; B-L, fixed 
transition 
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Figure 66: Skin friction envelope for the PAK B turbine stage - Re = 27, 000; k-e model 
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Figure 67: Skin friction envelope for the PAK B turbine stage - Re = 40, 000; B-L, turbulent 



Figure 68: Skin friction envelope for the PAK B turbine stage - Re = 40,000; B-L, floating 
transition 
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Figure 69: Skin friction envelope for the PAK B turbine stage - Re = 40, 000; B-L, fixed 
transition 



Figure 70: Skin friction envelope for the PAK B turbine stage - Re = 40, 000; k — e model 


85 





Figure 71: Skin friction envelope for the PAK B turbine stage - Re = 60, 000; B-L, turbulent 



Figure 72: Skin friction envelope for the PAK B turbine stage - Re = 60, 000; B-L, floating 
transition 
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Figure 73: Skin friction envelope for the PAK B turbine stage - Re = 60, 000; B-L, fixed 
transition 



Figure 74: Skin friction envelope for the PAK B turbine stage - Re = 60, 000; k - t model 
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Figure 75: Skin friction envelope for the PAK B turbine stage - Re = 80, 000; B-L, turbulent 



Figure 76: Skin friction envelope for the PAK B turbine stage - Re = 80,000; B-L, floating 
transition 
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Figure 77: Skin friction envelope for the PAK B turbine stage - Re = 80, 000; B-L, fixed 
transition 
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Figure 78: Skin friction envelope for the PAK B turbine stage - Re = 80, 000; k - t model 
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Figure 79: Skin friction envelope for the PAK B turbine stage - Re = 120, 000; B-L, turbulent 



Figure 80: Skin friction envelope for the PAK B turbine stage - Re = 120, 000; B-L, floating 
transition 
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Figure 81: Skin friction envelope for the PAK B turbine stage - Re — 120,000; B-L, fixed 
transition 



Figure 82: Skin friction envelope for the PAK B turbine stage - Re = 120, 000; k-e model 
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Figure 83: Skin friction envelope for the PAK B turbine stage - Re = 200, 000; B-L, turbulent 
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Figure 84: Skin friction envelope for the PAK B turbine stage - Re = 200, 000; k — e model 
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Figure 85: Momentum thickness envelope for the PAK B turbine stage - Re — 27,000; B-L, 
turbulent 



Figure 86: Momentum thickness envelope for the PAK B turbine stage - Re = 27, 000; B-L, 
floating transition 
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Figure 87: Momentum thickness envelope for the PAK B turbine stage - Re — 27,000; B-L, 
fixed transition 



Figure 88: Momentum thickness envelope for the PAK B turbine stage - Re = 27, 000; k — c 
model 
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Figure 89: Momentum thickness envelope for the PAK B turbine stage - Re = 40, 000; B-L, 
turbulent 



Figure 90: Momentum thickness envelope for the PAK B turbine stage - Re = 40,000; B-L, 
floating transition 
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Figure 91: Momentum envelope for the PAK B turbine stage - Re = 40,000; B-L, fixed 
transition 



Figure 92: Momentum thickness envelope for the PAK B turbine stage - Re = 40, 000; k — e 
model 
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Figure 93: Momentum thickness envelope for the PAK B turbine stage - Re = 60, 000; B-L, 
turbulent 



Figure 94: Momentum thickness envelope for the PAK B turbine stage - Re = 60, 000; B-L, 
floating transition 
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Figure 95: Momentum thickness envelope for the PAK B turbine stage - Re = 60,000; B-L, 
fixed transition 



Figure 96: Momentum thickness envelope for the PAK B turbine stage - Re = 60, 000; k — e 
model 
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Figure 97: Momentum thickness envelope for the PAK B turbine stage - Re = 80, 000; B-L, 
turbulent 



Figure 98: Momentum thickness envelope for the PAK B turbine stage - Re = 80, 000; B-L, 
floating transition 
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Figure 99: Momentum thickness envelope for the PAK B turbine stage - Re = 80, 000; B-L, 
fixed transition 



Figure 100: Momentum thickness envelope for the PAK B turbine stage - Re = 80, 000; k — t 
model 
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Figure 101: Momentum thickness envelope for the PAK B turbine stage -Re — 120,000; 
B-L, turbulent 



Figure 102: Momentum thickness envelope for the PAK B turbine stage - Re = 120, 000; 
B-L, floating transition 
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Figure 103: Momentum thickness envelope for the PAK B turbine stage - Re = 120, 000; 
B-L, fixed transition 



Figure 104: Momentum thickness envelope for the PAK B turbine stage - Re = 120, 000; 
k — e model 
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Figure 105: Momentum thickness envelope for the PAK B turbine stage - Re = 200, 000; 
B-L, turbulent 



Figure 106: Momentum thickness envelope for the PAK B turbine stage - Re = 200, 000; 
k — e model 
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Figure 107: Instantaneous entropy contours for PAK B turbine stage, Re = 40, 000, turbulent 
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Unsteady entropy contours for PAK B turbine stage, Re = 120, 000, turbulent 
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